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Abstract 

A new approach to the concept of particles and their production in quantum field theory is 
developed. A local operator describing the current of particle density is constructed for scalar 
and spinor fields in arbitrary gravitational and electromagnetic backgrounds. This enables 
one to describe particles in a local, general-covariant and gauge- invariant way. However, the 
current depends on the choice of a 2-point function. There is a choice that leads to the local 
non-conservation of the current in a gravitational or an electromagnetic background, which 
describes local particle production consistent with the usual global description based on the 
Bogoliubov transformation. The most natural choice based on the Green function calculated 
using the Schwinger-DeWitt method leads to the local conservation of the current, provided 
that interactions with quantum fields are absent. Interactions with quantum fields lead to 
the local non-conservation of the current which describes local particle production consistent 
with the usual global description based on the interaction picture. 



1 Introduction 



One of the main problems regarding quantum field theory in curved spacetime is how to introduce 
the concept of particles. The lack of a natural choice of the time coordinate implies that the 
generalization of the conventional definition of particles in Minkowski spacetime is not unique in 
general spacetime |2|. There are three main approaches to coping with this problem. Each of 
them has some advantages and disadvantages. Let us shortly describe them, emphasizing their 
disadvantages. 

The first approach is based on the theoretical point of view that the concept of particles has 
no fundamental meaning in field theory Q. Therefore, one should study only local, well-defined 
covariant operators, such as T^^ix), and try to express all observable quantities in terms of 
these. A disadvantage of such an approach is the fact that particles are what we observe in 
experiments. If we require that quantum field theory describes the observed objects, then it 
should describe particles. 

The second approach, complementary to the first one, resides on the physical point of view 
that particles are what particle detectors measure. Therefore, one avoids introducing a particle- 
number operator by studying the response of a model of a particle detector. These models 
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are usually of the Unruh-DeWitt type ^. However, such an approach is unable to answer 
the following questions. If all observables in quantum mechanics are represented by hermitian 
operators and do not require a model of a corresponding detector, then why is the particle 
number an exception? Do Unruh-DeWitt type detectors describe the essential properties of real 
particle detectors (such as a Wilson chamber or a Geiger-Miiller counter) in real experiments? 

The third approach, perhaps the most popular one, is somewhere between the first two ap- 
proaches. Although there may not exist a natural choice of the time coordinate in the whole 
spacetime, some regions of spacetime may have a natural choice of the time coordinate. Accord- 
ingly, one introduces different definitions of the particle-number operator for different regions. 
The two definitions are related by a non-trivial Bogoliubov transformation, which is interpreted 
as particle production |§, ^, |^] or as a variant of the Unruh effect |8|, |9|. In this way, one 
is able to talk about the average number of particles in some large region of spacetime. How- 
ever, the average number of particles should make sense in any region of spacetime, especially 
if the wavelength of the particles is much smaller than the region. The third approach, based 
on the Bogoliubov transformation, is not local and general-covariant, so it is unable to answer 
questions such as the following. How to calculate, at least in principle, the average number of 
particles in an arbitrary region of spacetime? Is the average number of particles in the union of 
two non-intersecting 3-dimensional regions that lie on the same 3-dimensional spacelike hyper- 
surface equal to the sum of the average number of particles in separate regions? If black holes 
radiate particles as predicted by Hawking 0, then what is the average particle distribution at, 
for example, 2 Schwarzschild radii away from the black-hole centre? 

There is a lot of similarity between the particle production by a background gravitational 



field and that by a background electromagnetic field |10| , |11| , 12 1. In particular, just as the 
former has problems with general covariance, the description of the later through the Bogoliubov 
transformation has problems with gauge invariance ||l^. This is related to the fact that gauge 
invariance is a local property, while the Bogoliubov-transformation method is based on a global 
definition of particles. 

Recently, a new approach to the concept of particles in quantum field theory has been 
proposed In this paper we further develop the ideas introduced in |14] and show that the 
new approach avoids the disadvantages of the other existing approaches discussed above. This 
new approach is inspired by some advantages of the first and the third approach and does not 
rely on any model of a particle detector. The main idea is the assumption that there should 
exist a well-defined local operator of particle-number density. There are two reasons for such 
an assumption. First, this allows a general-covariant and gauge-invariant formulation. Second, 
although the particle-number operator is a global operator in the conventional approach, this 
assumption is consistent with the fact that particles appear as local objects in all existing 
experiments. 

Our approach is based on a similarity between the number of particles and charge. For 
complex fields, the total number of particles is the sum of the number of particles and antipar- 
ticles, while the total charge is the difference of these two numbers. The concept of charge 
can be described in a local and covariant manner because there exists a local vector current 
of charge density. We find that a similar vector current exists for the number of particles as 
well. Nevertheless, it appears that this local current is not unique, but depends on the choice 
of a 2-point function. When a unique (or a preferred) vacuum exists, then the 2-point function 
is equal to the Wightman function and the current is conserved. When such a vacuum does 
not exist, then there is a choice of the 2-point function that leads to the local non-conservation 
of the current in a gravitational or electromagnetic background, which describes local particle 
production consistent with the usual global description based on the Bogoliubov transformation. 
Another choice, based on the Green function calculated using the Schwinger-DeWitt method 
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|T5| , p^ j, seems to be more natural. This choice leads to the local conservation of the current in 
an arbitrary gravitational or electromagnetic background, provided that other interactions are 
absent. 

Before studying the new local concept of particles in Sees. and [s], in Sec. § we clarify 
some physical aspects of the concept of particles in the conventional approach. In Sec. ^ we 
study in detail the case of hermitian scalar fields. In Sees. ^ and |5| we generalize the analysis to 
complex scalar fields and spin-1/2 fields, respectively. The conclusions and prospects for further 
investigations are given in Sec. ^. In the Appendix, we derive the particle-density operator in 
non-relativistic quantum field theory. 

2 Qualitative remarks on the concept of particles 

In quantum field theory, an n-particle state is defined as a normalized state of the form 

|n) = J d^h---d'knf{ki,...,K)al^---aljO). (1) 

The crucial question we attempt to answer in this section is why such formally defined states 
correspond to the states observed in typical experiments as n separated entities, i.e., "particles". 

We introduce the concept of a classical particle detector, such as a Wilson chamber or a 
Geiger-Miiller counter. We call such detectors classical, because, in order to understand how 
and why such detectors respond, quantum field theory is not essential. Classical detectors 
respond to states that correspond to the classical concept of particle, i.e., to states which are 
well localized in space. For example, if two localized particles are very near each other, which 
can be achieved by a suitable choice of the function /(ki,k2), then a classical detector will see 
this state as one particle. The response of a classical detector is also localized in time, in the 
sense that it can, with good precision, determine the instant of time at which the particle takes 
a certain position in space. In principle, the product of uncertainties AEAt may be arbitrarily 
small Because of these local features, one does not expect that the response of a classical 
detector depends on global aspects such as the whole detector's trajectory or the existence of 
an event horizon. Therefore, classical detectors are very different from particle detectors of the 
Unruh-DeWitt type. Because of locality, one expects that classical particle detectors detect 
objects that are covariant with respect to general coordinate transformations. Therefore, one 
needs a local and general-covariant notion of a particle in quantum field theory. 

As we mentioned above, a many-particle state can be detected as one particle if these lo- 
calized particles are very near each other. However, it seems that such a situation does not 
occur in practice. Why? Before answering this question, first consider how a one-particle state 
is detected. Assume that a particle detector is localized in space. The function /(ki) may 
correspond to a plane wave or to a state with two lumps, but such states cannot be detected 
by the localized detector. Such a state will not be detected until the function /(ki) collapses 
to a one-lump state. (We do not attempt to answer the question whether the wave function 
collapses spontaneously or the detector causes the collapse. In some interpretations of quantum 
mechanics there is a clear answer to this question, but we do not want to prefer any particular 
interpretation.) 

Now consider a two-particle state. Assume that, a priori, all functions /(ki,k2) are equally 
probable. (This probability does not reflect a stochastic nature of quantum states, but simply 
our ignorance of the details that govern the dynamics and the initial conditions.) However, 
the number of functions corresponding to a two-lump state is much larger than that for a one- 
lump state, so it is very improbable that two particles will form one lump. (Of course, the 
numbers of two-lump states and one-lump states are both infinite, but their ratio is not equal 
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to 1. To visualize this, note that if there are places where a lump can be placed, then (for a 
fixed number of particles) there are N different one-lump states and N'^ — N different two-lump 
states.) A four-lump state, for instance, is even more probable than a two-lump state, but if 
one particle is distributed in more than one lump, it will not be detected. 

In this way we have explained why an n-particle state behaves as an n-lump state in experi- 
ments; a particle splited in more than one lump will not be detected, while a smaller number of 
lumps can be realized in a smaller number of ways, so it has a smaller probability. Of course, if 
there are strong attractive forces among fields, there may exist a natural tendency for n-particle 
states to form a one- lump state (hadrons, a-particles), in which case it is more convenient 
to treat such states as one-particle states. Actually, it is often completely incorrect to treat 
such states as n-particle states, because non-perturbative effects of interactions may completely 
change the spectrum of states in a free theory, such as is the case for QCD. 

One may argue that classical detectors are not the best operational way to define a particle. 
One should rather study quantum detectors. (Even the response of a "classical" detector should 
be ultimately described by quantum field theory.) For example, an atom may be viewed as a 
quantum photon detector as it absorbs precisely one photon, not a half of it, nor two photons. 
It is often thought that this is an inherent quantum property. However, it is important to 
emphasize that there is not any deep, fundamental principle that forbids absorption or emission 
of a half of a particle or two particles. It is merely a consequence of a particular form of dynamics. 
For example, an electron in the atom absorbs one photon because the interaction Lagrangian 
-^int = ei^'JiiipA^ is linear in A^. Of course, higher-order corrections allow absorption of two 
photons (sum of their energies must be equal to the difference of energies of the atom levels) , but 
such processes are suppressed dynamically (small coupling constant) and kinematically (small 
probability of a one-lump two-photon state). However, is it, in principle, possible to absorb a half 
of a particle? One can exclude such a possibility by proposing that the Lagrangian C{(j), d^cj), . . .) 
of fundamental fields should be an analytic function around zero. For example, assuming that 
there are no derivative couplings, the analyticity implies that a local interaction Lagrangian has 
a form 

£int(0(x),---) = ECn(x)(</.(x))", (2) 
ra>0 

where n are integers and Cn{x) depend on some other fields. The interaction (§) implies that the 
number of absorbed or emitted particles must be an integer. If (j){x) is an effective, composite 
field, like (f){x) = xix)xix), then one can have a term proportional to \/0 = Xi i^i which case 
a "half of the 0-particle, i.e., one x-P^-rticle, can be absorbed or emitted. If one allows terms 
like for fundamental fields, then even a "half of a fundamental particle can be absorbed or 
emitted. (In this case, the concept of particle based on perturbative calculations is no longer a 
good concept. In particular, it is not clear, even algebraically, what VaJ\0) is.) 

3 Hermitian scalar field 

The hermitian scalar field is the simplest example of a quantum field. Therefore, we start the 
quantitative analysis by studying this simplest case in detail. The first two subsections do not 
contain new results, but they are written for the sake of completeness and easier understanding 
of our later results. 
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3.1 Canonical quantization 

Let g^u be a classical background metric, g the determinant of ^^j^, and R the curvature calcu- 
lated from g^^i^. The action of a hermitian scalar field (p is 

5 = 1 d^x\g\'/'^[g^'^id,<P){d,^)-m^^^-m% (3) 
where ^ is a coupling constant. Writing this as 

S = jd'x\g\'I^C, (4) 

the canonical-momentum vector is 

dC 

The correponding equation of motion is 

(V^a^ + m2+^i?)0 = O, (6) 

where is the covariant derivative. Let S be a spacelike Cauchy hypersurface with a unit 
vector normal to S. The canonical- momentum scalar is defined as 

TT = n'^TT/,. (7) 

The volume element on S is 

d^^" = d^x\g^^'>\^/^n^'. (8) 

The scalar product is defined as 



{<Pi,<P2) = i / d^f^r^d^cP^, (9) 

where ad^b = ad^h — {d^a)b. If cpi and 02 are solutions of @, then (^) does not depend on E 

We choose coordinates (t,x) such that t =constant on S. In these coordinates 

n'^ = (10) 
and the canonical commutation relations are 

[^{x),^{x')]^ = [7r(x),7r(x')]s =0, (11) 

[<A(x),7r(x')]s = |ff(=^)|-i/2i<53(x-x'). (12) 

The label S denotes that x and x' lie on S. Eq. ([l^ ) can be written in a manifestly covariant 
form as 

' dS'^[0(x), 5Xx')]x(x') = / d^'^^[(t>{x'),d^4>{x)]x{x') = ix{x), (13) 



where x(^0 is an arbitrary test function. 

For practical calculations, it is more convenient to introduce the quantity defined as 

= \gi3)\l/2^^l^ (14) 
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where the tilde on n'^ denotes that it is not a vector. It has the property V^n,/ = and has the 
exphcit form 

n'^ = (19^=^)^/^5^,0,0,0). (15) 
The volume element (^ can be written as 

dS^ = d=^xn^, (16) 

while (H) can be written as 

n^{x'mx),d'oct>{x')]j:=iS'{^-^'). (17) 
This can also be written in the canonical form as 

[<A(x),7f(a;')]s = ^<5'(x-x'), (18) 

where 

n=\g^^^\'/\ (19) 

is not a scalar. 

3.2 A choice of representation 

Let us choose a particular complete orthonormal set of solutions {fk{x)} of Eq. They satisfy 

{fk,fk') = -{fk^fk') = ^kk', 

{rkJk') = ifk,rk>) = o. (20) 

The field (j) can be expanded as 

^ix)=Y.^Mx) + alr,ix). (21) 



From (gD and we find 

afc = (/fc,c/>), al = -{ft,4>), (22) 
while from ( [20| ) and the canonical commutation relations it follows that 

[ak,ak'] = [4' 4'] = 0- (23) 

Therefore, we can interpret and a|. as lowering and raising operators, respectively. They 
induce the representation of the field algebra in the usual way. The vacuum |0) of the corre- 
sponding Hilbert space is defined by 

afcIO) = 0. (24) 
The operator of the total number of particles is 

N = J2alak. (25) 

k 

We also introduce the function W{x,x'), here defined as 

W{x,x') = Y,fk{x)mx'). (26) 
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(Later we also study different definitions oiW{x,x').) Using (^), we find that it is a Wightman 
function 

W{x,x') = {0\^{x)^{x')\0). (27) 



From ( P6D it follows that 

W*{x,x') =W{x',x). (28) 
From the fact that and satisfy Eq. (^) we find 

{V^'d^ + m^ + CR{x))W{x,x') = 0, (29) 

(V'^a' + + ^i2(x'))VF(3;, x') = 0. (30) 



From ( p7| ) and the canonical commutation relations we find that and are functions such 
that 

W{x,x')\j, = W{x',x)\j,, (31) 

ao5^Ty(x,x')|E = dod'oW{x',x)\^, (32) 

n°5^[VF(x,x') - VF(x',x)]s = i<53(x-x'). (33) 



3.3 The current of particle density 

The number operator given by ( p5D is a global quantity. However, as originally found in |14| j, a 
new way of looking into the concept of particles emerges when ( p2|) is put into (|25| ) , and (^)and 
( p6| ) are used. This leads to 

N = j^dS^ dS'^VF(x, x') a^ai, 0(x)0(x'). (34) 

By interchanging the names of the coordinates x and x' and the names of the indices ^ and 
this can be written as a sum of two equal terms 

N = - f dJ:'' f dJ:"'W{x,x')df,dl(t>{x)(t>{x') 
2 Js Jt. 

+ \ I d-E^ I dT.'^W{x\x)d^l4>{x')'t){x). (35) 

2 

Using also (p8[), we see that ( |35| ) can be written in the final form as 

N = J^d^'^j^ix), (36) 
where the vector jf^ix) is the hermitian operator 

j^ix)= I dT.'^\{W{x,x')d^^U{x)<l){x') + h.c.}. (37) 
Js 2 

Obviously, the vector j^{x) should be interpreted as the local current of particle density. The 
representation of by (^) and ( |37| ) has three advantages with respect to the conventional 
representation (p5|). First, it avoids the use of the operators and a\ related to a particular 
choice of the modes fk{x)- Second, it is manifestly covariant. Third, the local current j^{x) 
allows to view the concept of particles in a local manner, which the conventional representation 
does not allow. 
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If we put ( piD and (Pq) in (p7|) and then use @, (pO|), (|2^), and the antisymmetry of the 
operators^, we find 

j^l = ^Yl fk d^i fk'a\ak' ■ (38) 

k,k' 

From this form we see that is automaticaUy normally ordered and has the property 

JmIO) = 0. (39) 

This is actually not surprising because we have started from the operator (25) which is also 
normally ordered. 

Using (^) and (p9|), we find that the current ( |37| ) possesses the property 

V^^j.^ix) = 0. (40) 

This covariant conservation law means that the background gravitational field does not produce 
particles, provided that a unique (or a preferred) vacuum defined by (^) exists. (This local 
covariant conservation implies also the global conservation because it provides that ( |36|) does 
not depend on time. The extra terms in (^0|) that originate from the fact that ^ are 
compensated by the extra terms in (|3^ ) that originate from the fact that dS/^ / (d^x, 0, 0, 0).) 
The choice of the vacuum is related to the choice of the function W{x, x') (see (|2^ and (p^)). 

Note that although i^(x) is a local operator, some non-local features of the particle concept 
still remain, because (|37|) involves an integration over S on which x lies. Since (t){x') satisfies (^) 
and W{x,x') satisfies (30), this integral does not depend on S. However, it does depend on the 



choice of W{x,x'). Note also that the separation between x and x' in (|37|) is spacelike, which 
softens the non-local features because W{x,x') decreases rapidly with spacelike separation. As 
can be explicitly seen with the usual plane- wave modes in Minkowski spacetime (see Sec. |3.4| ), 
W{x, x') is negligible when the spacelike separation is much larger than the Compton wavelength 

3.4 The case of Minkowski spacetime 

In Minkowski spacetime with the metric g^y = diag(l, —1, —1, —1), it is most natural to choose 
the modes fk as the usual plane- wave modes /k, where k is the 3- momentum. Using the discrete 
notation for the summation over k, these modes are 

where = Vk^ + m^, V is the 3- volume and = (L<Jif,k). Eq. (Bsl) becomes 



The physical content is more obvious from the space-integrated operator 

J^.^ j d^xj^ix) = Y.Zr "k«k' (43) 



k 



where the tilde on denotes that it does not transform as a 4- vector. The quantity k^/uj]^ is 
the 3- velocity u', so we find 

J\nq) = VUglUq), (44) 
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where \nq) = (n^!) ^/^(aj)""'|0) is the state with Uq particles with the momentum q. Eqs. (|4 
and ( |4^ ) support the interpretation of as the particle current. 

-ik{x—x') 



The two-point function (2f:) is 



k 2u;k 

In the infinite- volume limit, we perform the replacement 

k 

SO (H) takes the usual form 



(27r)3' 



(45) 



(46) 



l^(x,x') 



-ik{x—x') 



(47) 



The properties of this function, often denoted by A(3; — x') |T^, are well known. In particular, 
for t = t' and |x — x'| = r, it behaves asymptotically as 



a/ TTIT 

A{r) oc ^ e""^'' for rm > 1. 



(48) 



For m = 0, the function W{x,x') is 
A(x - x')\m=o 



1 

1^ 



1 



7 ^ + i7r5((x -x')^) 

(x — x'j^ 



(49) 



The function ( |47| ) is also equal to G'^{x,x'), where iG'^{x,x') is equal to the Green function 



G{x,x')- J ^2^^^ ^2_^2 ' 
calculated with the appropriately chosen contour of integration over ko . 
3.5 The question of positivity of the local number of particles 



(50) 



Eq. ( |36D can be written as 



where 



N = J d^a; 1^(3) 1 1/2^(2;), 



n = 



(51) 



(52) 



Since the total number of particles ( psD is a non-negative operator (in the sense that the expected 
value {ip\N\tp) is non-negative for any state IV')), one could expect that the local density of 
particles n{x) is also a non-negative operator. For example, in the Minkowski case, for \q) = 
OqlO), we find {q\n\q) = {q\jo\q) = V~^i which is a positive constant. However, for = 
2^"'^^^(ki) + \q2)), where qi / 92, we find 



(V'IJoIV') = 77 



1 + ^ cos(gi - q2)x 



(53) 



which is negative for some values of x, provided that / a;q2. Therefore, the particle density 
n(x) is not a non-negative operator. Is this in contradiction with the experimental fact that the 
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number of particles in a small volume cr C S cannot be negative? The answer is no! When the 
position of a particle is measured with some accuracy such that it is (almost) certain that the 
particle is localized inside o", then the state {ip) of the particle is a localized state, such that the 
corresponding wave function vanishes (or is negligible) outside of a. Consequently, the quantity 
{'ip\n{x)\ip) vanishes (or is negligible) outside of a. Therefore, the total number of particles in 
this state is (approximately) equal to 



(fx \g 



(3) 1 1/2 



(V'ln(x)lV'). 



(54) 



Since the total number of particles cannot be negative, it follows that the number of particles 
inside a is non-negative for such a localized state. 

Let us illustrate the heuristic arguments above by a concrete example. We study the case of 
Minkowski spacetime with infinite volume V , so the field expansion (21) becomes 



4){x) 



d^A;[ak/k(x)+aj^/i^(x)], 



where 



-ikx 



(55) 
(56) 



Let the state of the system be the one-particle state 

|y)^c-i/V(0,y)|0). (57) 
The constant c is the norm determined by the requirement that (y|y) = 1, so it is equal to 

c = A+(0) = A~(0), (58) 

where 



A±(x) 



(27r)32LJk 



(59) 



The wave function corresponding to the state (|t)Y|) is 



%by{x) = c{0\(l){x)\y) = cc-^/2^+(t, X - y), (60) 

where c is a constant chosen such that the wave function is normalized in an appropriate way. 
This wave function is related to a Lorentz- invariant notion of particle detection |20|. Using the 
fact that W{x,x') = A^{x — x'), a straightforward calculation gives 



(y|jo(2;)|y) 



where 



A+(t,x-y)(^-(t,x-y) 
A+(0) 



c.c. 



(2^)2 



In particular, at t = 



(y|jo(0,x)|y) =<53(x-y), 



(61) 



(62) 



(63) 



which implies that the state (57) has a non-negative and strictly localized particle density at 
t = 0. Note that the wave function (60) is not strictly localized at t = 0, but decreases rapidly 
with |x - y| (see (||) and (||)). 
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Note also that in a typical physical situation in which particles are detected (see Sec. |2|), the 
state {tp) corresponds to particles, where each particle is localized inside one of lumps. In 
such a situation, ( p4[ ) represents the number of particles in one of the lumps, so this number is 
equal to one. 

We also stress that negative particle densities do not appear in non-relativistic quantum field 
theory (see Appendix A). Heuristically, this can be seen from the frequency-dependent factor 
that multiplies the cosine function in (^). Since in the non-relativistic limit cjg ~ m, this factor 
becomes equal to 1, leading to a non- negative density. 

Note finally that, in general, the particle density (that may be negative) is not the same 
thing as the probability density (that cannot be negative) of finding the particle at a position 
X. The probability density in a state is equal to |(x|^)p, where |x) is an eigenstate of the 
position operator X with the eigenvalue x. For the relativistic case, the problem of finding the 
operator X and the states |x) is not yet satisfactorily solved. In the nonrelativistic limit, the 



states |x) are equal to the states |x) defined by (57). 



3.6 Particle production by non-gravitational interactions 

Let us now study the case in which a non-gravitational interaction represented by £int = — Wint 
is also present. In this case, the equation of motion is 

{Vd^ + + ^R)(j) = J, (64) 

where J{x) is a local operator containing cj) and/or other dynamical quantum fields. Since it 
describes the interaction, it does not contain terms linear in quantum fields. (For example, when 
£int = — A(/>^/4, then J(x) = —X(j)'^{x)). We propose that even in this general case the particle 
current is given by ( p7| ) 

j^ix) = J^di:"'^{Wix,x')d^dl^{x)^{x') + h.c.}, (65) 

where W{x, x') is the same function as before, satisfying the "free" equations 

{V^'d^. + m^ + iR{x))W{x,x') = (66) 
{V'^d'^ + + iR{x'))W{x, x') = 0, (67) 

and having the expansion 

W{x,x')=Y^h{x)ft{x'). (68) 

k 

As we show below, such an ansatz leads to particle production consistent with the conventional 
approach to particle production caused by a non-gravitational interaction i2int- 
The field (j) in (65) satisfies (p3). Therefore, using (|§), we find 



(x)= f dJ:"'-{W{x,x')dlJ{x)(t>ix')+]i.c.}. (69) 



Note that only x (not x') appears as the argument of J on the right-hand side of ([69|), implying 
that J plays a strictly local role in particle production. 



Let us now show that our covariant description (69) of particle production is consistent with 
the conventional approach to particle production caused by a non-gravitational interaction. Let 
T,{t) denote some foliation of spacetime into Cauchy spacelike hyper surf aces. The total mean 
number of particles at the time t in a state is 

Mit) = imitm, (70) 
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where 



N{t) = / dS'^j^(t,x). (71) 

JS(t) 



Equation ( [70D is written in the Heisenberg picture. However, matrix elements do not depend 
on picture. We introduce the interaction picture, where the interaction Hamiltonian is the part 
of the Hamiltonian that generates the right-hand side of (|6^). The state 1-0) transforms to a 
time-dependent state 

IV'int(t)) = Uintitm. (72) 

Here the unitary operator Uintit) satisfies 

i^Uintit) = Hintit)Uintit), (73) 

where 

^int(t) = / d^x\g\'/^nUx)- (74) 

Note that Eqs. ([7^ ) and ( |7^ are not manifestly covariant. The covariance of the time-evolution 
law is discussed in Sec. |3.9| . In the interaction picture, the field (j) satisfies the free equation (^), 
so the expansion (^l|) can be used. This, together with (ffll), (|65|) , and (|68|), implies that N{t) 
transforms to 

k 

Therefore, (^) can be written as 

M{t) = (V'mt(t)|iVi„t|V'int(i)) = I](0mt(t)l4afcl^mt(i)>, (76) 

k 

which is the usual formula that describes particle production caused by a quantum non-gravitational 
interaction. This is most easily seen when t — > oo and the solution of ( |73D is chosen such that 
|V'int(i — *■ —oo)) = \il>), in which case C/int is the ^-matrix operator. 



3.7 Particle production by the gravitational background 

So far, our discussion has been based on a particular choice of the modes {fki^)} as preferred 
modes. In practice, they are chosen such that fk are positive frequency solutions of (^). However, 
when metric is time dependent, then such preferred modes do not exist. In the conventional 
approach to the concept of particles, this is related to particle production by the gravitational 
background [§, ^. In this case, one introduces a new set of functions ui{x) for each time t, such 
that ui{x) are positive- frequency modes at that time. This means that the modes ui possess an 
extra time dependence, i.e., they become functions of the form ui{x;t). These functions do not 
satisfy (^). However, the functions ui{x\t) satisfy (^), provided that r is kept fixed when the 
derivative acts on ui. To describe the local particle production, we take 

W{x,x') = Y.ui{x;t)ul{x'-t'), (77) 

instead of (^). Since ui{x;t) do not satisfy (|6|), the function ( [77| ) does not satisfy (p9|). Instead, 
we have 

{Vdf, + + ^R{x))W{x, x') = -K{x, x') + 0. (78) 
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Using (0) and (^) in ([37|), we find a relation similar to (f 



V%(x)= / d^"'l{K{x,x')dl^{x)^{x')+h.c.}. (79) 
This local description of particle production is consistent with the usual global description 



based on the Bogoliubov transformation. This is because and ( |37D with (|77| ) and (21) lead 
to 

Nit) = J2Alit)Mt), (80) 

where 

Mt)=Y,aUt)ak-P!kit)al (81) 

k 

mk{t) = ifk,ui), Pikit) = -Ului). (82) 

The time dependence of the Bogoliubov coefficients aik(t) and Pik{t) is related to the extra 
time dependence of the modes ui{x;t). If we assume that the change of the average number of 
particles is slow, i.e., that 

dtMt) ^ 0, (83) 

which occurs when 

dtUi{T,^;t)\r=t^O, (84) 

then the Bogoliubov coefficients (^) are equal to the usual Bogoliubov coefficients. This approx- 
imation is nothing else but the adiabatic approximation, which is a usual part of the convential 
description of particle production |Q . 

3.8 The natural choice of W{x, x') 

In general, there is no universal natural choice for the modes ui{x;t). In particular, in a given 
spacetime, the choice of the natural modes ui{x;t) may depend on the observer. If different 



observers (that use different coordinates) use different modes for the choice of (77), then the 
coordinate transformation alone does not describe how the particle current is seen by different 
observers. In this sense, the particle current is not really covariant. There are also other problems 
related to the case in which different observers use different modes [pl| . 

Since covariance was our original aim, it is desirable to find a universal natural choice of 
W{x,x'), such that, in Minkowski spacetime, it reduces to the usual plane- wave expansion (|47|). 
Such a choice exists. This is based on the fact that the Feynman Green function Gf{x, x') can be 
calculated using the Schwinger-DeWitt method and that the knowledge of the Feynman Green 



function automatically yields the knowledge of other Green functions 16|, such as G'^{x,x'). 
The Schwinger-DeWitt method does not require a choice of a particular set of modes. Therefore, 
we propose 

W{x,x') = G+{x,x'), (85) 
where G"*" is calculated from Gp given by 

/•oo 

GF{x,x ) = i[g{x)g{x')]-^l^ ds(x, s|x', 0). (86) 

JO 

The function (x, six', 0) is determined by the function cj(x, x') and its derivatives, where (j{x, x') 
is one half of the square of the geodesic distance between x and x' [^] . Therefore, the Green 
function Gi? is unique, provided that a geodesic connecting x and x' is chosen. When x and 
x' are sufficiently close to one another, then there is only one such geodesic. In this case, the 
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adiabatic expansion |2|, ^] oi Gp can be used. If Riemann normal coordinates are spanned 
around x', then all Green functions can be written as |2|, ^] 

G{x,x') = \g{x)\-'/'J ^e-'=(— (87) 



where 



g{k;x') = Qn{k]x') 

n=0 



and Qn{k;x') have n derivatives of the metric at x' . The two lowest functions of the adiabatic 
expansion (^) are the same as in Minkowski spacetime, i.e., Qq = [k"^ — m?)~^, Qi = 0. De- 
pending on the choice of the integration contour over fco, from (|87|) one can calculate Gp, G"*", 
or any other Green function. In particular, G^ satisfies 

{V^'d^ + n? + ^R{x))G+{x,x') = (89) 

{V^d'^ + + iR{x'))G^{x, x') = 0. (90) 

Since G~^{x, x') decreases rapidly with spacelike separation, the contributions to (^) from large 
spacelike separations may be negligible. Therefore, for practical calculations, it may be sufficient 
to calculate G^{x,x') only for small spacelike separations, in which case the first few terms of 
the adiabatic expansion are sufficient. 

When x and x' are not close to one another, then there may be more than one geodesic 
connecting x and x' . In this case, (x,s|x',0) depends on the choice of the geodesic, the point x 
is not contained in the subspace covered by the Riemann normal coordinates spanned around x' 
and the adiabatic expansion does not work. Nevertheless, (x,s|x',0) can be computed for any 
choice of the geodesic connecting x and x' . The resulting Green function G^ satisfies (|89| ) and 
(pO|). To define the exact unique particle current, we need a natural generalization of G^ to the 
case with more than one geodesic connecting x and x' . The most natural choice is the 2-point 
function defined as the average over all geodesies connecting x and x' . Assuming that there 
are N such geodesies, this average is (for any Green function G) 

N 

G(x,x') = iV-^^G(x,x';aJ, (91) 

a=l 

where G{x,x' ;aa) is the Green function G{x,x') calculated with respect to the geodesic Ga- 
in the case with a continuous set of geodesies connecting x and x', the generalization of 
( |9lD is more complicated. At the point x' we span the flat Minkowski tangent space. Then we 
perform the Wick rotation on the tangent space which transforms it to a Euclidean tangent space. 
Each geodesic emanating from x' is uniquely determined by the space angle O = ((^, "i^i, i92) in 
the Euclidean tangent space because the space angle determines the direction of a geodesic 
emanating from x' . Therefore, the measure on the set of all geodesies emanating from x' is (see, 
for example, p3|]) 

/ dVL^> = [ dip r sm^id^i f sm^'&2d^2, (92) 
J Jo Jo Jo 

where geodesies emanating in two opposite directions are counted as two different geodesies. Of 
course, not all geodesies emanating from x' cross the point x, but some of them do. The Green 
function corresponding to a geodesic determined by its space angle ilx' at x' and crossing the 
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point X is denoted by G{x,x'; ilx'). Therefore, we introduce the averaged Green function 



dVtx'G{x, x' ] ^x') 
Gx'{x,x') = ^ , (93) 



where /' denotes integration only over those space angles for which the corresponding geodesic 
crosses the point x. However, a geodesic connecting x and x' can also be specified by its angle 
in the Euclidean tangent space spanned at x, resulting in the averaged Green function 



dj^xG{x ^ X 5 

Gx{x, x') = J, . (94) 



In general, ( p3| ) and (p^ are not equal. To treat the points x and x' in a symmetric way, we 
introduce the symmetrized Green function 

G{x,x') = ]^[Gx{x, x') + Gx' {x,x')]. (95) 

The symmetrized Green function satisfies (p8|), i.e. 

G*{x,x')=G{x',x), (96) 

provided that G is a Green function such that G*(x,x') = Gx{x',x) and G*,{x,x') = Gxi{x',x). 
Now the most natural choice for the particle current in arbitrary curved background is 

j^{x) = |^dS'^i{G+(x,x')5,S,0(x)(^(x') + h.c.}. (97) 



Since the function G^{x,x') satisfies ( [8^ ) and (pO|), it follows that the particle current (97) does 
not depend on the choice of E and is conserved 

V^i^(x) = 0, (98) 

provided that the field satisfies (P). Note that the corresponding definition of particles does 
not always correspond to the quantities detected by "particle detectors" of the Unruh-DeWitt 
type 0, ^. Instead, the number of particles resulting from ( p7| ) is determined by a well-defined 
hermitian operator that does not require a model of a particle detector, just as is the case for 
all other observables in quantum mechanics. Moreover, since the observable (/'(x) in (|97| ) does 
not require a choice of representation of field algebra, the definition of particles based on (5+ 
does not require the choice of representation either. This allows us to treat the particles in the 



framework of algebraic quantum field theory in curved background |24]. Furthermore, as noted 
by Unruh [Q, only one definition of particles can correspond to the real world, in the sense 
that their stress-energy contributes to the gravitational field. Since the definition of particles 
based on G^ is universal, unique, and really covariant, it might be that these are the particles 
that correspond to the real world. Similarly, from the discussion in Sec. |2|, it seems that these 
particles might correspond to the objects detected by real detectors (such as a Wilson chamber 
or a Geiger-Miiller counter) in real experiments. If such an interpretation of these particles 
is correct, then classical gravitational backgrounds do not produce real particles, in agreement 



with some other results |T^, |2^, |2^, Note that this does not necessarily imply that 
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black holes do not radiate. It is possible that quantum gravity provides a correct mechanism for 
black hole radiation. Whatever the mechanism of the particle creation near the horizon might 
be, one should expect an approximately thermal distribution of escaped particles as seen by a 
distant observer. The thermal distribution is a consequence of the exponential red shift and 
can be understood even by classical physics |2^, without any assumption on the physical 
mechanism that causes particle creation near the horizon. 



3.9 Covariance of the time evolution 

As we have seen, when particle production occurs, then the function W{x,x') or the field 4>{x') 
does not satisfy the "free" equation of motion ( |30| ) or @, respectively. Consequently, the 
current j^{x) depends not only on the point x, but also on the choice of S on which x lies. 
This suggests that there exists a preferred foliation of spacetime. Such a result should not be 
surprising because many other aspects of quantum theory seem to require the existence of a 
preferred time coordinate. However, it may be surprising that, in our analysis, the need for 
a preferred time coordinate has emerged from the requirement of general covariance (of the 
concept od particles). Below we argue that the principle of general covariance (together with 
some other natural principles) may help us to determine this preferred time coordinate. 
The time evolution law (|7^ ) written in the Heisenberg picture is 

i^U{t) = H{t)U{t), (99) 

where (see (^)) 

H{t) = [ (fx\g\^/^n{x), (100) 

n = T^, and 

Similarly, in the Schrodinger picture 

ij^^[(t),t) = M[<^,t), (102) 
where ^[(p^t) is a functional with respect to (/'(x) and a function with respect to t, while the 



Hamiltonian H is given by (100) in which ■7r(x) = —i5 / 5(f){x.) (see (|l8|)). This time-evolution law 
is not manifestly covariant. This is because t and H are not scalars. This law would become 
covariant if dt H could be replaced with dr H, where 

dT = n.dx'' = -^dx^", (103) 

H= f d^>'n''T.,= f d'x\g^^^\'/^^^ (104) 
Jn 500 

are scalars. It is easy to see that dtH = dTH for any Tj} if and only if 

QQi = and diQQQ = 0. (105) 



In this case, dr = \/goo{t)dt, so the metric takes the form 

^ _ g^,^^^ x)dx'dx^. (106) 
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The metric can always be written in this form. The coordinates for which the metric takes the 



form ( |106| ) are known as Gaussian coordinates |31]. 



The results above can be reinterpreted in the following way. The covariant time evolution 
determined by (|103| ) and ( |104D depends on the foliation of spacetime into Cauchy spacelike 
hypersurfaces. The choice of foliation corresponds to the choice of coordinates for which t = 
constant on S. Therefore, we can determine the time evolution uniquelly by giving a preferred 
status to a set of coordinates. The Gaussian coordinates are special because only in these 
coordinates the time-evolution law takes the canonical form (]99|), required by the usual rules of 
quantum mechanics. This shows that the requirement of covariance of the time-evolution law 
leads to Gaussian coordinates as preferred coordinates. 

However, for a given spacetime, the Gaussian coordinates are not unique. One needs an 
additional criterion for choosing them. For example, a reasonable additional criterion is the 
requirement that there should be no coordinate singularity in these coordinates, but this re- 
quirement still does not make the Gaussian coordinates unique. An interesting attempt to 



choose the Gaussian coordinates in a unique way is given in |32, However, in the ap- 



proach of [32, 33], the coordinate singularities such as horizons are not avoided. In the case 
of a Schwarzschild black hole, a particularly interesting choice of the Gaussian coordinates is 
the Lemaitre coordinates, in which there is no coordinate singularity. In these coordinates, the 
Schwarzschild black hole is time dependent, the Hawking radiation looks thermal to a distant 
outside observer, but the absence of the horizon implies a unitary time evolution which resolves 
the black- hole information paradox |34]. 

4 Complex scalar field 

To explore the similarity between the current of particles and that of charge, we generalize 
some results for hermitian scalar fields to the case of complex scalar fields. We do not discuss 
those aspects the generalization of which is obvious. We also study the scalar QED, i.e., the 
interaction of complex scalar fields with classical and quantum electromagnetic fields. 

4.1 Particle and charge currents 

The complex field (p can be written in terms of hermitian fields (pi and 02 as 

(pi + i(p2 



V2 ■ 



(107) 



The independent fields (pi and (p2 satisfy the same canonical commutation relations as the 
hermitian scalar field of Sec. ^, with an additional property [(pi{x),(p2{x')] = 0. Therefore, the 
canonical commutation relations can be written as 

[(P\x), = [do(p\x), d'o(P{x')h = 0, (108) 

n°(x')[(/>^(x), 9^0(x')]e = iSH^ - x'). (109) 
The complex field (p satisfies (|6|) and can be expanded as 

k 

(P\x) = Y,alfl{x) + hMx). (110) 
k 
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The operators (a^) destroy (create) particles, while bk destroy (create) antiparticles. 
They satisfy the corresponding algebra of lowering and raising operators and can be expressed 
as 

ak = {fkA), h{ = -{fl,4)), 
ai = -{n,^^), bk = {fk,4>^)- (111) 
We introduce two global quantities: 

N^^^ =J2alak±blbk. (112) 

k 

Here A^^"*") is the total number of particles, while A^^"-* is the total charge. In a way similar to 
the case of hermitian field, we find the covariant expression 

ivW= l^dS^jW(x), (113) 

where 

j!^H^) = |^dS'-^{t^(x,x05^i[c/>1-(:E),/.(:E') ±</.(x),/.t(x')] +h.c.} (114) 

and W{x,x') is given by (p6|). Obviously, j^^^ should be interpreted as a current of particles, 
while ^ should be interpreted as a current of charge. These two currents can also be written 



as 

,{±) _ ,(P) ^ , 

where 



^=jP±jif\ (115) 



jPi^) = lj^"'l{Wix,x')d^dlcp\x)^{x') + h.c.}, (116) 
A^Hx) = I dT.'^]-{W{x,x')d^dlct>{x)ck\x')+h.c.} (117) 



/e 2 

are the currents of particles and antiparticles, respectively. It is also instructive to write the 
currents in terms of hermitian fields and 02 . The current of particles is a sum of two currents 
of the form (37); 

= |^dI]"^^{t^(x,x')5.S[<^i(x)</.i(x') +</>2(x)02(x')] +h.c.}, (118) 

while ili ^ takes a form from which it is manifest that the current of charge does not exist in the 
case of only one hermitian field; 

= |^dS'^^{VF(x,x')5X«['^i(^)'/'2(a:') -02(x)0i(x')] +h.c.}. (119) 

There is one important difference between the currents j/j^"* and \ Similarly to the 
particle current (|37|), the particle current j\j^\x) possesses the non-local features related to 
the integration over S on which x lies. On the other hand, the apparent non-local features of 

\x) really do not exist, because, by using the canonical commutation relations ( [1081 )- ( p^9D 

+) 



and Eqs. (^)-(|3^), the integration over dTi'^ can be performed. This cannot be done for j/t 
because, owing to the different sign, certain terms that cancel in ^ do not cancel in j^^'^ ■ 
Applying the commutation relations ( 108 )-( |l09|) to ^ given by (|114|) such that (f)^ always 



18 



comes to the left and cf) to the right, a straightforward calculation that exploits Eqs. (|3l|)-(33) 
yields 

j\:^\x) = i4)\x)d^(t){x) + [ <m"'w{x,x')d^lw{x',x), (120) 



so all non-local features are contained in the second term that does not depend on (f). Similarly, 
applying the commutation relations such that (f)^ always comes to the right and cj) to the left, 
^ can be written as 

= - / dT."'W{x\x)diSlW{x,x'). (121) 



From the antisymmetry of the operators and 5^, it follows that the integrals appearing in 
(|12C| ) and (|121| ) are equal. Therefore, by summing ( |120D and ( |12lD and dividing the sum by 2, 
we find a purely local expression 

= ^[0t(x)5.(/.(x) (122) 

This is the usual form of the charge current, often derived as a Noether current resulting from 
the global U(l) invariance of the Lagrangian. 

Using (|11[1|), one can show that (|122|) is normally ordered, i.e., that 



(0|j(-^|0) = 0. (123) 
Alternatively, one can prove ( |123| ) by using ( |110| ) , (|^) , and (^0|) to show that 

(o|i0^(x)al(/)(x)|o) = - / (n:"'w{x,x')dj^lw{x',x), (124) 

which implies that (|120| ) can be written as 

Using (H), (|TlO|) , (|0|), and [0^,6^'] = 0, the currents (|ll|), ([Tl^ ), and ( p7|) can be written 
in a form similar to (|38|). We find 



J'u'=^T. f*k df. A'4«fc' + (126) 



k,k' 



j!^^ = ^j:rAfk4h'-jr, (127) 

k,k' 

ji^^=iJ2fAh'i4^k' + blbk'), (128) 

k,k' 

Jt^ = iT^fAfM^k' - blbk>) + 2jf^, (129) 



k,k' 



where the mixed hermitian current mixes (in products) the particle operators Ofc, with 
the antiparticle operators bk, b^: 

= I EifAh'ak'bk + fAfkAbl)- (130) 

k,k' 



19 



Now we see that 

Ji+^|0)=0, (131) 

while the currents ( |126| ), ( |127| ), and ( |129| ) have not this property. In the case of Minkowski 
spacetime, Eq. ( |43| ) generalizes to 

j{±) ^ J d'xj^tKx) = Y: ^(4«k ± ^l&k). (132) 
From (^) and (^) it fohows that the currents (|114| ) are conserved: 



V^i^) = 0. (133) 

(+) 



When a non-gravitational interaction is described by (^4|) , then ' satisfies an equation similar 
to (H), i.e. 

V^j(+)(a;)= / (n:'''-{W{x,x')dl[j\x)<P{x') + J{x)4>\x')]+h.c.}. (134) 
Js 2 

On the other hand, ^ is conserved again, provided that the Lagrangian possesses a global 
U(l) symmetry. The conservation of ^ is related to the fact that the second term in ( |120| ) 
does not depend on the interaction described by J. 

When a natural choice of the modes ft does not exist, then one can redefine the particle 
current ( |114D by replacing (p6|) with (77) or G^{x,x') from Sec. 



4.2 The generalization to the electromagnetic background 

Let us generalize the definition of currents j^^^ to the case in which a classical electromagnetic 
background A^{x) is also present. Under a gauge transformation 



the scalar field transforms as 



4(x) = A^{x) + e-i9^A(x), (135) 
4,\x) = (t>{x)e-'^^'-'l (A'^x) = </.t(x)e*^(^). (136) 



The covariant derivatives 

D^ = V^ + ieA^, Dl = V^-ieA^ (137) 

are covariant with respect to both the gauge transformations and the general-coordinate trans- 
formations. The quantities D^j_(j) and D*^(f)* transform as and cj)* in (|136| ), respectively. The 
fields (f) and (j)* satisfy the equations of motion 

{D^'D^ + m2 + ^R)(f) = 0, 

+ m2 + ^i?)</)t = 0. (138) 

The solution can be expanded as 

<^(^) = H akfk{x) + blglix), 

k 

(/.t(x) = ^ 4a*(x) + bkgkix). (139) 

A; 
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When Afj_ = 0, then fk = Qk- In general, fk and g^. transform as (j) under a gauge transformation, 
while and gk transform as cj)^ . Introducing the "antisymmetric" covariant derivative!)^ defined 

by 

aDf,b = aD^b- {Dla)b, (140) 
we define two different scalar products: 



(a, 6) =i f d^''a*D^b, 



{a,b)^ = i J d^''a*D*b. (141) 

These scalar products are gauge invariant, provided that a and b transform in the same way 
under the gauge transformation, i.e., either as (j) or as (j)^ . When a commutes with b, then 

aD*^b = -bDf,a. (142) 
Consequently, in the commuting case, the two scalar products are related as 

{a,b), = -{b*,a*). (143) 

The modes fk and gk satisfy 

{fk,fk') = -{fk^fk')* = ^kk', 
{9k, 9k')* = -{ghdk') = ^kk', 

{fk,9l,) = {9k,rk')*=0- (144) 

Therefore, the operators ak and bk satisfy the usual algebra of lowering and raising operators 
and can be written as 

ak = {fk,4>), bl = -{gl,(j)), 
4 = -(/fe'<^^)*' bk = {9kA^)*- (145) 

The vacuum |0) is defined by 

ak\0) = bk\0) = 0. (146) 
We introduce two 2-point functions: 

w^pHx,x') = Y.fkix)mx'), 

k 

W(^\x,x') = J2gk{x)gUx'). (147) 

k 

From ( 139 ) and the fact that the operators Uk and aj, are related to particles while bk and 6|. are 
related to antiparticles, we see that W^^'^ is related to the propagator of particles, while W^^^ is 
related to the propagator of antiparticles. In general, W^^^ ^ W^^^ , which corresponds to the 
fact that particles and antiparticles propagate in a different way, provided that a background 
electromagnetic field exists. When = 0, then W^^^ = W^^^ = W, where W is the 2-point 
function used in Sec. 4.1. This means that particles and antiparticles propagate in the same 
way when a background electromagnetic field does not exist. When F^j, = dfj,A,^ — d^A^j^ = 
but Afj_ = e^^d^X / 0, then 

W^P\x,x') = e"^[^(^)-^(^')lH^(a;,a;'), 

W^^\x,x') = e^[^(^)-^(^')lVF(x,x'), (148) 
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so the two 2-point functions differ only by a phase, which should not have any physical effects. 
The functions W^^^ and W^^^ satisfy 

{D'^D^ + m^+ iR{x))W^^\x, x') = 0, 
{D'*^'D'* + m^ + iR{x'))W^P\x, x') = 0, 
{D*i'Dl + m? + iR{x))W^^\x, x') = 0, 

{D'l'D'^ + m2 + (x, x') = 0. (149) 

Therefore, VF*-^-* and W^^^ transform to each other when transforms to —A^^ which corre- 
sponds to the transformation F^^ — > —F^^. From (139) it follows that VF*^^) and W^^'^ are equal 
to the Wightman functions 

W^P\x,x') = (O|0(2;)0t(x')|O), 

W^^\x,x') = (O|(/.t(x)0(x')|O), (150) 

while from (147) we find 

W^P>{x,x') = W^^\x',x), W^(^)*(x,x') = W^^\x',x), (151) 

Under gauge transformations, W^^\x,x') transforms as 4'{x)<j)^ {x'), while W^'^^x^x') trans- 
forms as (t)\x)(t){x'). From the canonical commutation relations 

[0t(x),0(x')]E = m\x),D'^^{x')h = 0, (152) 
h\x')[(t>^{x),D'^4>{x')]^=i5\^-^'), (153) 
we find that and gt are functions such that 

[W^P\x,x')-W'^^\x\x)]^ = Q, 
DoD'q* [W^(^) {x, x') - {x', x)]e = 0, 
h^{x')D'o*[W^^\x,x') - W^^\x' ,x)]j: = i6^{^ - ji.'). (154) 

Now we introduce the global quantities: 

N^^^ =J2alak±blbk. (155) 
k 

Writing this in the covariant form (|113|), we find 



jW(x) = J^d^'^^{[W^PHx,x')D;Dl^Hx)cP{x') 

±W'^'^\x,x')Df,D'; (l){x)(l)\x')] + h.c.}. (156) 

Although VF(^), W^^\ (/>, and (/>t depend on gauge, the currents ( |156| ) do not depend on gauge. 
Therefore, our definition of particles and that of charge are gauge invariant. Eq. (156) is the 
gauge-invariant generalization of (114). 

In a way similar to the case = 0, it can be shown that ^ can be written in a purely 
local form 

= -icPD;<p^ + {0\i<pD;^^\0). (157) 
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This is the usual form of the charge current, often derived from the gauge invariance of the 
Lagrangian. 



From ( |149| ) and (138) it follows that both currents are conserved: 

V^jj^^ = 0. (158) 

Of course, the conservation of the charge current jji ^ was expected because it is well known 
that (157) is conserved. However, the conservation of jj^^ is less trivial because it implies 
that the background electromagnetic field described by does not produce particles (particle- 
antiparticle pairs), provided that a unique (or a preferred) vacuum defined by (146) exists. 

In general, a natural definition of the vacuum ( p,46| ) does not exist because the natural choice 
of the modes fk and gk does not exist. In this case, one can introduce different modes that are 
natural at different times, similarly as in Sec. 3.7. In a time-dependent background the 
natural modes at a given time can be chosen such that they are positive-frequency modes at 
that time ||T^. Alternatively, in a time-independent but space-dependent background A^^, one 
can introduce "in" and "out" modes, the naturalness of which is related to the sign of the wave 
vector k in asymptotic regions [BHI. In this way, one can describe the pair creation by the 



electromagnetic background in a local way, similarly as in Sec. |3.7| . As in the case of particle 
creation by the gravitational background studied in Sec. it is easy to see that this local 
description of pair creation is equivalent to the usual global description based on the Bogoliubov 
transformation. 

However, the sign of the frequency and that of the wave vector are not gauge-invariant 
concepts. In some cases, such as F^u = or a time-independent background electric field, 
some gauges may seem "natural", but a general recipe for choosing the gauge does not exist. 
Therefore, by using the Bogoliubov-transformation method, it is not possible, even in principle, 
to calculate the distribution of produced pairs in an arbitrary electromagnetic background. 

The following remark also illustrates why the Bogoliubov-transformation method is prob- 
lematic. For F^i, = 0, the "natural" gauge is = 0, which defines the "natural" vacuum. 
However, by using a different gauge, one can obtain that this vacuum is a many-particle state. 
In particular, an appropriate gauge leads to a thermal distribution of particles in the "natural" 
vacuum, leading to an electromagnetic analog of the Rindler quantization and the Unruh effect. 
On the other hand, physical quantities should not depend on gauge, which again raises doubts 
about the appropriateness of the Bogoliubov-transformation method. 

Now we see that just as the particle current based on (^^ is not really general-covariant, a 

similar modification of ( |156| ) is not really gauge- invariant. To define particles in a unique way, 

one needs a universal natural choice of the functions W^^^ and W^^^ using a method that does 

not require a particular choice of the modes fk and g^- Again, such a choice exists and is a 

I I (p\ 

generalization of the choice in Sec. 3.8. By introducing the Feynman Green functions Gp and 

Gj^^ calculated by using Schwinger's method jl^, 10, |l^, we propose 

W^^\x,x') = G(^)+(x,x'), W^'^\x,x') = G(^)+(x,rc'). (159) 

(Note that the Green functions are calculated for Minkowski spacetime in the references cited 
above, but the formalism can be generalized to curved spacetime as well.) Here G^^^~^ and 
satisfy (|l49|), i.e. 

{D^'D^ + m2 + ^fl(x))G(^)+ (x, x') = 0, 
{D'*^'D'* + m2 + iR{x'))G^P^+{x, x') = 0, 
(D*^'Dl + + iR{x))G^^^+ (x, x') = 0, 

{D'^'D'^ + + ^i?(x'))G(^)+ (x, x') = 0, (160) 
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and are calculated from the corresponding Feynman Green functions that satisfy similar equa- 
tions, but with the right-hand sides proportional to S^{x — x'). Note that in the literature 
1 15, 10, 11, p^], usually only one Feynman Green function corresponding to clf^ is considered 



(P) 

and is called simply Gp. However, there are actually two Feynman Green functions Gp and 
G^p^ that satisfy similar equations, but with replaced with —A^. 



From (|159|) , (|16C| ), and (|138|) it follows that the currents ( |156|) are conserved. In particular, 
this implies that the most natural definition of the particle current in an electromagnetic back- 
ground leads to the result that classical electromagnetic backgrounds do not produce particles. 



This is in agreement with some other results 1 26, 27, 25]. It is interesting to note that a different 



gauge-invariant formalism based on the same Green function Gp as above leads to a result that a 



classical background static electric field produces particle-antiparticle pairs |15]. However, when 
the formalism of ||l^ is naturally modified so that the integration contour over a scalar variable 
remains on the real axis (see (^6|)), then this formalism does not lead to pair production. As 
argued in |25], such a modified formalism might be the correct one. 



4.3 The generalization to the quantum electromagnetic interaction 

As we have seen, the result that the interaction with a classical electromagnetic background does 
not lead to a change of the number of particles is closely related to gauge invariance. However, 
it is known from many experiments that the electromagnetic interaction may cause a change 
of the number of charged particles. The best known example is the annihilation of a particle- 
antiparticle pair to a pair of photons. Such processes are successfully described by Feynman 
diagrams that result from a theory in which the electromagnetic field is also quantized. Physical 
arguments explaining why only a quantum electromagnetic interaction can consistently describe 



a change of the number of particles are given in |28|. Here we give a formal explanation of this. 
Let A^ be the total operator of the electromagnetic field. It can be written as 

A^ = A^ + A^, (161) 

where is the operator such that 

{0\A^\0) = 0, (162) 

from which it follows that 

(0|i^|0) = A^. (163) 

When Afj^ does not interact with other fields, then it can be expanded in the usual way in terms 
of creation and annihilation operators. A gauge transformation 

A'^ = A^ + e-'d^X, (164) 

where A(x) is a c-number function, can be written as 

A' — A 

A'^ = A^ + e-'d^X. (165) 

Therefore, the operator A^ is invariant with respect to the classical gauge transformations ( |165| ). 
The generalization of ( |137| ) is 

D^ = Vf, + ieAf, = Df, + ieA^. (166) 
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The currents j/i take the same form (|156| ) with the same covariant derivative and with the 
same functions W^^^ and VF^"^) satisfying Eqs. ( |149| ). However, in these currents, the field (f) 
satisfies 

{b^'D^ + m^ + ^R)(t) = 0, (167) 



and similarly for Eq. (|167| ) can also be written in a form similar to (|64|): 

(D'^L)^ + + ^R)<j) = J, (168) 

where 

J = [e^A'^^^ - ie(V^A'')](/> - lieAf'D^cj). (169) 



We see that J transforms under the gauge transformations ( |165| ) in the same way as (/>, so ( |168| ) 
is covariant. Similarly, the currents jji^^ given by ( |156| ) are invariant with respect to gauge 
transformations ( |165| ). The local particle production is described by the gauge-invariant gener- 
alization of (|134| ). In a similar way as before, it can be shown that this local particle production 
is consistent with the particle production described in the usual way using the interaction pic- 
ture, in which the interaction Hamiltonian is the part of the Hamiltonian that generates the 
right-hand side of (|168| ). 

The conservation of charge is a little bit tricky. One can introduce two charge currents: 

q^ = i4>^D^(t>, q^ = i(P^D^(P. (170) 

Again, it can be shown that 

jl"^=9M-(0kM|0)=:9M- (171) 



However, this current is not conserved. From the equation of motion (167) it follows that the 
current is conserved, while is not: 

V^g^ = 0, 

V^g^ = 2e[{V ^,A^')cp^ + A^'V ^{(t>^)]. (172) 



Note that the right-hand sides in (|172|) do not depend on the background A^, so the following 
peculiarities appear even when A^ = 0. The two currents are related as 

= - 2e^^0V. (173) 

Which charge current is more physical, q^ or q^l It depends on how the charge is measured. 
If the charge is measured so that the corresponding electromagnetic effects determined by the 
total electromagnetic field F^^ = d^Ay — d^Afj, are measured, then the current is the physical 
one, because the Maxwell equations read 

V''F^, = -e:q^: . (174) 

However, this does not mean that the current q^ is unphysical. Each hermitian operator, con- 
served or not, is an observable. In order to understand the physical meaning of q^, in the rest 
of this subsection we work in the interaction picture, in which the field (j) has the expansion 
(|13S| ) and A^ has a similar expansion in terms of creation and annihilation operators. If one 
simply counts the number of particles and the number of antiparticles, then the current q^ is 
the physical one because the corresponding charge 

Q= j d-^^:q^:= N^-^ = ^ a^a^ - foj^ftfc (175) 
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is equal to the difference of these two numbers. On the other hand, the conserved charge 

Q = J dW ■.q^:= N^-^ -2e J d^'^Af, ■.^^4>: (176) 

depends not only on the number of charged particles, but also on the number of photons. 

An interesting question is the following. Can one observe, by counting the number of scalar 
particles and antiparticles, that the charge Q is not conserved? Since the extra term in ( |173| ) 
is linear in the expected value of this operator vanishes for any state in which the number 
of photons is definite (i.e., the state is not a superposition of states with different numbers of 
photons). Therefore, the expected values of the charges Q and Q are equal for such states, so 
the charge Q is conserved in this case. Note also that the number of any kind of particles is 
definite in any particular Feynman diagram, so the charge Q is also conserved on the level of 
Feynman diagrams. The conservation of Q on the level of Feynman diagrams can also be seen 
directly from the fact that, in the interaction Hamiltonian, the fields (f> and (j)^ always appear in 

products of the form and (j)^ D^cj). However, if the initial state \i) is a state with a definite 
number of all kinds of particles, then the final state determined by the S'-matrix operator takes 
the form 

\f) = S\i) = Y^C,f^U). (177) 

n 

Here the states | fn) have a definite number of all kinds of particles, but the total state | /) has not. 
(The coefficients Cj/^ = {fn\f) = {fn\S\i) are the amplitudes determined by the corresponding 
Feynman diagrams.) Therefore, it is possible, in principle, to find the system in a state in which 
the final charge Q is not equal to the initial charge Q. However, it seems very unlikely that 
this will happen in practice because Q is conserved in any experiment in which the state |/) 
collapses to a state with a definite number of particles, or the decoherence destroys all mixing 
terms so that the observable degrees of freedom can be described by a density matrix that is 
not a pure state, but a classical mixture of the states with different numbers of particles. 

Note, finally, that analogous problems with the conservation of charge do not exist for spin- 
1/2 fields. In the Standard Model, fundamental charged scalar fields do not exist. There is the 
complex Higgs field, but, after the spontaneous symmetry breaking, it reduces to an uncharged 
hermitian Higgs field. Therefore, the problems with the conservation of charge for complex 
scalar fields may be phenomenologically irrelevant. 

5 Spinor field 

For simplicity, we study the spinor fields in Minkowski spacetime. However, all equations can be 
generalized to an arbitrary spacetime by replacing the derivatives with the general-covariant 
derivatives and by replacing certain constant Dirac matrices with the x-dependent Dirac matrices 
7^(3;) = e^(x)7a, where e^(x) is the vierbein 0. When appropriate, we shortly discuss the effects 
of this generalization. 

5.1 Particle and charge currents 

The Dirac matrices satisfy the anticommutation relations 

{7;„7.} = 2r?^.- (178) 
We choose their representation such that |jl^ 

7^ = 7o7m7o- (179) 
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The massive spin-1/2 field tplx) satisfies the Dirac equation 

(i ^ - = 0, +m) = 0, (180) 
where = j^d^ and '0 = ■0^70- The field can be expanded as 

'^i^) = ^kUk{x) + d\vk{x), 

k 

'^(^) = X! A'^kix) + dkVk{x), (181) 



where the summation over k includes the summation over spin indices as well. The scalar 
product is 

(V,x) = / c^S'^V^T^X- (182) 

JT. 

The spinors Uk and Vk satisfy 

{uk,Uk') = ivk,Vk') = Skk', 
{uk,Vk') = {vk,Uk') = 0. (183) 

Therefore, 

bk = {uk,ip), dl = {vk,ip), 

bl = {ij,Uk), dk = {il^.Vk). (184) 
The canonical anticommutation relations can be written as 

{'^a{x),'4'h{x')]T. = {{'4>{x)lQ)a, {Hx')lo)b]T. = 0, 

{Va(x),^c(x')}(70)c6|E = <5afe5^(x-x'). (185) 



From ( |185| ) it follows that the operators (|184D satisfy the usual algebra of fermion creation and 
annihilation operators: 

{hk,hl,] = {dk,dl,] = 6kk', (186) 
while other anticommutators vanish. We introduce the particle and antiparticle 2-point functions 

S^^\x,x') = ^Uk{x)Uk{x'), 
k 

S^^\x,x') = Y.Vk{x')vk{x). (187) 



They satisfy 



{i0-m)S^^\x,x') = 0, S^^\x,x'){i p' +m) = 0, 

5(^)(x,a;')(i0 +m) = 0, (i - m)S^^\x,x') = ^ (188) 



and have the property 



(5(^)(x,x'))a6 = (0|V'a(x)V5fe(x')|0), 

(S(^)(x,x'))a6 = mh{x)^aix'm. (189) 
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We also introduce the 2-point functions s\!^){x,x') and s\^\x,x') defined as 

S\f,^^\x,x')=^,S^P'^\x,x')^,. (f90) 

Using (|17|), (|l79l) and dTs^) we find 

^oS^P'^'^\x,x')^o = S^P'^\x\x\ (191) 

while from ([T78|) , ([l79|) , ( |l85|) , and (jl8|) it follows that 

[S^P){x,x')- S^^\x\x)]^=-to5^{^-^'), 
[S\l!^\x,x') - S\:^\x' ,x)]y. = 7^7o7i.5^(x - x'), 

[S%\x,x') - S^^\x\x)]^ = 7,<53(x - x'). (192) 

The global quantities 

^^^^=E^K±44 (193) 

k 

can be written in the covariant form as 

ivW = /^dS'^J?), (194) 

where 

j^^\x) = j^d^'^\mx)S\!^\x,x')^{x') 

±V^(x)4^)^(x, x'){f^{x')] + h.c.} (195) 

and the superscript T denotes the transpose. From ( |180 ) and ( |18§| ) it follows that both currents 
are conserved: 

S'^ii^^ = 0. (196) 

By an appropriate generalization to a curved background, it is straightforward to see that the 
currents j^^^ are conserved again. Also, one can introduce different natural spinors at different 
times, which redefines the 2-point functions (187) and leads to a local description of particle pro- 
duction by the gravitational background consistent with the Bogoliubov transformation method. 
The most natural choice of the 2-point functions is 

S^P\x,x') = {ij/ + m)G+{x,x'), 

5(^)(x,x') = G+(x,x')(iy -m), (197) 

where G^{x,x') is the 2-point function for scalar fields with ^ = 1/4 discussed in Sec. 
This natural choice also leads to the conservation of the particle current. 
Using ( |185| ) and ( [L92[ ) we find the standard form of the charge current 



i(-)(x) = 4>{x)^^i^{x) - Tr5(^)(x,x)7^, (198) 
which, after using ( |181| ) and ( |187D , can be written as 

= i>{x)lM^) - (0|V5(x)7^7/'(a^)|0). (199) 
The currents j^^^ can also be written as 

i^)=jr±i^^. (200) 
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In terms of the creation and annihilation operators, the currents can be written as 



jP = J2^knf.Ukbl,bk + jf^ (201) 

k,k' 



i^^ = T.^m'dldk-jr, (202) 

k,k' 

= "^k'li^Ukhl/hk + Vk7i,Vk'dl,dk, (203) 



k,k' 



(-) 



where 



^(nfc'7^tifc6^,6fc - Vk'jfj,Vk'dl,dk) + 2j™'', (204) 

k,k' 

3^"" = \ ^{vkltiUk'dkbk' + UkniiVkb\,d\). (205) 

k,k' 



5.2 The generalization to the electromagnetic interaction 



The expressions for the currents in Sec. 5.1 are already gauge invariant. This makes the gen- 
eralization to the case of electromagnetic interaction easier than for the scalar fields. This is 
related to the fact that the Dirac equation is a first-order equation, so the expressions for the 
currents contain 7^ instead of d^. Consequently, in the currents, one does not need to replace 
the derivatives with the gauge-covariant derivatives. Therefore, the currents j'^^^ are given by 
the expression ( [L95| ). However, the fields satisfy 

{ip - m)ij = 0, ^{i p* +m) = 0, (206) 
while the 2-point functions satisfy 

{ip -m)S^^\x,x') = 0, S^^\x,x'){ip'* +m) = 0, 

S^^\x,x'){ip* +m) = {), (ip'-m)S(^)(x,x') = 0. (207) 



The charge current can be written in the form ( |199| ) again and is conserved. On the other hand, 
the particle current is not conserved: 

a^j(+)(x) = |^dS'-i{ieA^(x)[V;(x)4f)(x,x')V(x') 

-V^(x)5(^)^(x, x')V5^(x')] + h.c.}. (208) 

The background does not appear in (|208D , which means that the background electromagnetic 
field does not produce particles. Again, it may be interpreted as a consequence of the existence 
of the preferred spinors Vk and u/c, or as a consequence of the existence of the natural 2-point 
functions calculated by using Schwinger's method The particle production by the 

electromagnetic background can be described by introducing different natural spinors at different 
times, which redefines the 2-point functions and leads to a local description of particle production 
consistent with the Bogoliubov transformation method. The particle production described by 
(|20j) is equivalent to the usual description based on the interaction picture. 
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6 Conclusion and outlook 



In this paper, the currents of particles and charge for scalar and spinor fields in gravitational and 
electromagnetic backgrounds have been constructed. The currents are covariant with respect to 
general coordinate transformations and invariant with respect to gauge transformations. The 
currents of charge constructed using our method are equal to the usual charge currents. However, 
previous definitions of particles in quantum field theory were not general-covariant and gauge- 
invariant. This is because previous definitions of particles were not local, i.e., the operator of 
the local particle density was not known, except for the non-relativistic case. 

For a given background, the particle current is not unique, but depends on the choice of 
a 2-point function. Different choices correspond to different definitions of particles. There are 
3 types of this choice. The first type is based on the choice of a particular set of complete 
orthonormal modes that satisfy the field equations of motion. The resulting particle currents 
give a local description of the particle content in various inequivalent representations of field 
algebra. The second type is based on choosing different natural modes at different times. This 
leads to a local description of particle production by classical backgrounds. When the total 
number of particles on a Cauchy hypersurface is calculated, then the first and second types of 
choice lead to two types of conventional global concepts of particles. However, the third type, 
based on the 2-point function calculated using Schwinger's method, is novel even on the global 
level. This type seems to be the most natural one. It is tempting to interpret these particles as 
real particles. If this interpretation is correct, then classical backgrounds do not produce real 
particles. There are also other indications that classical backgrounds might not produce particles 



|13| , 25, 26, 27, Even if these particles do not correspond to real physical particles in general, 
it is interesting to ask about the physical meaning of this hermitian observable that corresponds 
to physical particles at least in Minkowski spacetime and zero electromagnetic background. 

As by-products, relevant even without our novel description of particles, the following results 
have also been obtained. A unique Green function G{x, x') has been constructed for the case 
in which x and x' can be connected by more than one geodesic. It has been found that the 
requirement of covariance of the time evolution leads to a preferred role of Gaussian coordinates. 
The peculiarities related to the concept of charge and its conservation in scalar QED have been 
clarified. 

In this paper, the phenomenological applications of the developed formalism have not been 
studied, but we hope that our results will motivate further investigations. We suggest several 
possible applications. 

First, one can adopt the variant of the formalism according to which background gravitational 
and electromagnetic fields do produce particles. Contrary to the previous methods, our method 
enables one to calculate where particles are created. This question is particularly interesting in 
the case of particle production by black holes. Our formalism is able to test the usual conjecture 
that particles are created near the horizon. 

Second, one can adopt a variant of the formalism according to which background gravitational 
and electromagnetic fields do not produce particles. It seems very likely that this is the correct 
interpretation because it results from the definition of particles that seems to be the most natural. 
In order to calculate particle currents or the total number of particles on a spacelike Cauchy 
hypersurface, one has to calculate 2-point functions explicitly by using the Schwinger-DeWitt 
method, which is a non-trivial task in general. 

Third, one can study possible physical implications of the fact that the local particle den- 
sity may be negative. It would be interesting to study more thoroughly how the appearance of 
negative particle densities depends on the state, entanglement and background interactions, as 
well as how negative densities disappear by the wave-function collapse or decoherence. For ex- 
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ample, the appearance of negative particle densities could be related to the existence of classical 



"barriers" such as a horizon or an electric potential higher than 2m. Eq. (38) suggests that neg- 
ative particle densities may be related to positive-norm solutions that contain both the positive 
and the negative frequencies, but this depends on which of the 3 types of 2-point functions is 
adopted. 

In any case, we believe that our results could influence further research that could result in a 
significant change of the understanding of the concept of particles in quantum field theory. Our 
local approach changes the notion of a relativistic particle even for a free particle in Minkowski 
spacetime and zero electromagnetic background. 
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A Particle density in non-relativistic quantum field theory 

In non-relativistic quantum field theory, the operator of particle density is a well-known quantity 
p7| , This is a non-negative operator 

n{x) = ^'■^(x)^'(x), (209) 

where x = (x, t) and ^ satisfies the Schrodinger equation 

-^V^^ + U^ = t^, (210) 
2m ot 



with U{x) being a non-relativistic potential. In this appendix, we rederive ( P09D using the 
method employed in this paper. 

Eq. (|210| ) can be derived from the Lagrangian density 

C = i^^- (V^'1')(V*) - ^7^%. (211) 

Therefore, the canonical momentum is tt = i"^^ and the canonical equal-time commutation 
relation between ^ and vr can be written as 

[*(x,t),^t(x',t)]=53(^_^/), (212) 

The scalar product {ct)i,4'2) = / d^x(j)l(j)2 does not depend on time, provided that (pi and (j)2 
satisfy (|21^). The fields ^ and *t have the expansion 

f (x) = J2 akUkix), ^^(x) = alulix). (213) 

k k 

The modes Uk are orthogonal: 

{uk,Uk') = iul,ul,) = 6kk'- (214) 

From ( |213| ) and ( pl4D we find 

ak = iuk,^), 4 = (4,*t). (215) 
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From the equal-time commutation relations it follows that and o|. satisfy the usual algebra 



of lowering and raising operators jS^] , so the operator 

N = J2aW (216) 
k 

is the operator of the number of particles. From ( |213| ) we see that 

{0\¥{x)^{x')\0) = 0, (0|^'(x)^^(x')|0) = W{x,x'), (217) 

where 

W{x,x') = Y.Uk{x)uUx'). (218) 

k 



From ( pl2| ), ( 217 ), and (|21§| ) we see that the modes Uk are complete, i.e., that 

W^((x,t),(x',t)) = 53(x-x'). (219) 
Using ( |215|) and ( plSj) , we can write (|216|) as 

N = J (fxn{x), (220) 

where 

n{x) = J d^x'-{W{x,x')'^\x)'i'{x') + h.c.}. (221) 

In ( ^21]) , t = t', so (221) and ( ^19| ) lead to (|209D . Therefore, contrary to the relativistic case, in 
the non-relativistic case the particle density is non-negative and can be written in a purely local 
form. 

For comparison with the relativistic case, let us study the free case [7 = 0. The modes Uk 
can be chosen to be the plane-wave modes uj^, which, in a finite volume V, are 

e-«{^k*-kx) 

txk(x,t) = -= , (222) 



where = k^/2m. In particular, taking the state to be = 2 ^^^(|qi) + |q2))) which is the 
non-relativistic analog of the state \ip) in Eq. (|53|), we find 



(^InlV') = V-^{1 + cos[(£;qi - Eq,)t - (qi - q2)x]}. (223) 

Contrary to the relativistic particle density (|5^), the non-relativistic particle density (|223| ) is 
non-negative. 
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